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ABSTRACT 

Although they are but a small fraction of the mass ejected in core-collapse supernovae, neutrino- 
driven winds (NDWs) from nascent proto-neutron stars (PNSs) have the potential to contribute 
significantly to supernova nucleosynthesis. In previous works, the NDW has been implicated as a 
possible source of r-process and light p-process isotopes. In this paper we present time-dependent 
hydrodynamic calculations of nucleosynthesis in the NDW which include accurate weak interaction 
physics coupled to a full nuclear reaction network. Using two published models of PNS neutrino 
luminosities, we predict the contribution of the NDW to the integrated nucleosynthetic yield of the 
entire supernova. For the neutrino luminosity histories considered, no true r-process occurs in the 
most basic scenario. The wind driven from an older IAMq model for a PNS is moderately neutron- 
rich at late times however, and produces ®^Rb, *^Sr, ^^Y, and ^"Zr in near solar proportions relative 
to oxygen. The wind from a more recently studied 1.27Mq PNS is proton-rich throughout its entire 
evolution and does not contribute significantly to the abundance of any element. It thus seems very 
unlikely that the simplest model of the NDW can produce the r-process. At most, it contributes to 
the production of the N = 50 closed shell elements and some light p-nuclei. In doing so, it may have 
left a distinctive signature on the abundances in metal poor stars, but the results are sensitive to both 
uncertain models for the explosion and the masses of the neutron stars involved. 
Subject headings: nuclear reactions, nucleosynthesis, abundances, supernovae — stars: neutron 



1. INTRODUCTION 

The site where r-process nuclei above A=90 have 
been synthesized remains a major u nsolved prob- 
lem in nucleosynthesis theory (e.g., lArnould et al.l 
|2007[ ). Histo rically, many possibilities have been 
proposed (see IMeveH |1994|) . but today, there are 
two principal contenders - neutron star mergers 
(jLattimer et al.l 119771: iFre iburghaus ct al. 1999) and 
the NDW fWooslcv ct al.' 1994:' Oian & Wooslcv 1996 



Hoffman et al. 1997; Otsuki ct al. 2000; Thompson ct al.l 
20011) . Observations of ultra-metal-poor stars suggest 



that many r-process isotopes were already quite abun- 
dant at early times in the galaxv dCowan et al.l 119951 : 
iSneden et all 119961 : iFrebel et all 120071 ). suggesting both 
a primary origin for the r-process and an association 
with massive stars. NDWs would have accompanied the 
first supernovae that made neutron stars and, depending 
upon what is assumed about their birth rate and orbital 
parameters, the first merging neutron stars could also 
have occurred quite early. 

Both the merging neutron star model and the NDW 
have problems though. In the simplest version of galac- 
tic chemical evolution, merging neutron stars might be 
capable of providing the necessary integrated yield of the 
r-process in the sun, but they make it too rarely in large 
doses and possibly too late to be consistent with obser- 
vations (Argast et al. 2004). On the other hand, making 
the r-process in NDWs requires higher entropies, shorter 
time-scales, or lower electron mole numbers, Yg, than 
have been demonstrated in any realistic, modern model 
for a supernova explosion (though see iBurrows et al.l 
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I2OOI . 

Previous papers and models for nucleosynthesis in the 
NDW have focused on the production of nuclei heav- 
ier than iron u sing cit her greatly si mplified dynamics 
( Beun et"all[20 b8: Farou oTet al.ll2009D or nuclear physics 
jPian & Wooslcv 1996; Otsuki ct al. 2000:'Arconcs et al.l 
I2OO7.: .Fischer ct al. 2009; Hucdcpohl ct al. 2009). Post 
processing nuclear network calculations have been per- 
formed using th ermal histories from accurate models 
of the dynamics ([Hoffman et al.l Il997t I Thompson et al.l 
IMntlWaMlol 120061) . but the calculations sampled only 
a limited set of trajectories in the ejecta. No one has 
yet combined the complete synthesis of a realistic NDW 
with that of the rest of the supernova. 

To address this situation, and to develop a frame- 
work for testing the nucleosynthesis of future explo- 
sion models, we have calculated nucleosynthesis using 
neutrino luminosity histories take n from two PNS cal- 
culations found in th e literature (jWooslev et al.l 119941 : 
iHuedepohl eTaIll2009l ). This was done using a modified 
version of the implicit one-dimensional hydrodynamics 
code Kepler, which includes an adaptive nuclear network 
of arbitrary size. This network allows for the production 
of both r-process nuclei during neutron-rich phases of the 
wind and production of light p-elements during proton- 
rich phases. Since the results of wind nucleosynthesis 
depend sensitiv ely on the neutrino luminosities and in - 
teraction rates (|Oian fc WooslevllT996l : iHorowitzl I200I . 
we have included accurate neutrino interaction rates that 
contain both general relativistic and weak magnetism 
corrections. 

The synthesis of all nuclei from carbon through lead 
is integrated over the history of the NDW and combined 
with the yield from the rest of the supernova, and the 
result is compared with a solar distribution. If a nucleus 
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produced in the NDW is greatly overproduced relative 
to the yields of abundant elements in the rest of the su- 
pernova, there is a problem. If it is greatly underpro- 
duced, its synthesis in the NDW is unimportant. If it is 
co-produced, the NDW may be responsible for the galac- 
tic inventory of this element. An important outcome of 
this study are the yields expected from a "plain vanilla" 
model for the NDW. Are there any elements that are ro- 
bustly produced and thus might be used as diagnostics 
of the wind in an early generation of stars? 

In iJJl we discuss the general physics of neutrino driven 
winds and analytically delineate the regions in neutrino 
temperature space were different modes of nucleosynthe- 
sis occur. We then discuss our numerical model in Sj3l 
In f|4l the results of the time dependent models are pre- 
sented. We conclude with a discussion of how the NDW 
might affect galactic chemical evolution and consider if 
this allows the strontium abundance in low metallicity 
halo stars to be used as a tracer of supernova fallback at 
low metallicity. Additionally, we investigate if observed 
abundances in SN 1987A can put constraints on late 
time neutrino luminosities from PNSs. Finally, we dis- 
cuss some possible modifications of the basic model that 
might improve the r-process production. These ideas will 
be explored more thoroughly in a subsequent paper. 

2. GENERAL CONCEPTS AND RELEVANT 
PHYSICS 

After collapse and bounce in a core collapse super- 
nova, a condition of near hydrostatic equilibrium exists in 
the vicinity of the neutrinospheres. The temperature of 
the outer layers is changing on a time scale determined 
by the Kelvin-Helmholtz time of the PNS, t kh ~ 10s 
(|Burrows k Lattimeri IToM iPons et al.lll999f ). This is 
much longer than the dynamical time scale of the PNS 
envelope, so the hydrostatic part of the envelope is in 
an approximate steady state. The neutrino heating rate, 
which is determined by the neutrino luminosities from 
the neutrino sphere, must then balance the local neu- 
trino cooling rate. Heating and cooling are dominated 
by the charged current processes (fe, + n) ^ (e~ + p) 
and {Pe+p) ^ {e++n) |Qian fc Wooslevlll996l ). Equat- 
ing these rates, while neglecting the neutron-proton mass 
difference and weak magnetism corrections and assuming 
the geometry can be approximated as close to plane- 
parallel gives the temperature structure of the neu- 
tron star atmosphere as a function of radius, Tatm ~ 
IMMeY R~l'^Ll[^^el'lj^y (y^/y)^^^, where L^,5i and 
ii^,MeV are the electron neutrino luminosity and aver- 
age neutrino energy at the neutrino sphere in units of 
Iq5i gj.gg g-i g^]-^j MeV, r espectively. The gravitational 
redshift factor is y — -^/l — 2GMns /rc^ which, when 
evaluated at the neutrino sphere, R^, is y^. Notice that 
the only dependence on radius is carried in the redshift 
factor, so that the atmosphere is close to isothermal. 

At the radius, Tc, where the pressure in the enve- 
lope becomes radiation do minated, the mater i al be - 
comes unstable to outflow (Salpet er fc Shapiro! 119811 ). 
Since the neutrino luminosity is significantly lower than 
the neutrino Eddington luminosity , a thermally driven 
wind results (jDuncan et al.l I1986D . The density at 
which this wind begins can be found approximately by 
equating the radiation pressure to the baryonic pres- 
sure. This results in a critical density, pc ~ 8.3 x 



10'^ gem '■^ R^^QL]J'^-^€i,MeV {yv/y) , at which significant 
outflow begins and the kinetic equilibrium of weak inter- 
actions ceases to hold. Under these conditions, nuclear 
statistical equilibrium is maintained on a time scale much 
shorter than the dynamical time scale and, for these tem- 
peratures and densities, there will be no bound nuclei 
present. Since the electron fraction is set by kinetic equi- 
librium, the composition of the wind does not depend on 
any previous nuclear processing, so any nucleosynthesis 
from the wind will be primary. 

Assuming that most neutrino heating occurs near r^, 
the entropy is constant once the temperature cools to 
the nucleon recombination temperature, kT k, 0.5 MeV. 
Therefore, the final nuclear abundances in the wind de- 
pend mainly on the wind entropy, electron fraction, and 
the dynamical ti mescale at the radius w h ere alpha com- 
bination occurs ()Qian fc WooslevI 119961 : iHoffman et"all 
Il997| ). To determine the contribution of the wind to the 
nucleosynthesis of the entire supernova, the mass loss 
rate must also be known. Estimates for these quantities 
are given in the Appendix along with a discussion of the 
effect of general relativistic corrections. 

Integrating the mass loss rate (equation lAlip for a typ- 
ical neutrino luminosity history implies that the wind 
will eject approximately 10~^ Mq of material. This in 
turn means that for the wind to contribute to the inte- 
grated yields of the supernova for a particular isotope, 
that isotope needs to overproduced relative to its solar 
mass fraction by a factor of at least 10^ in the wind, as- 
suming the rest of the supernova ejects 10 Mq and has 
over production factors of its most abundant metals of 
order 10. 

At early t imes in the wind, the PNS is losing lepton 
number (e.g. lBurrows &: Lattimeilll986( ). Neutrino inter- 
actions in the wind then tend to increase the lepton num- 
ber of the wind, and, to maintain charge neutrality, cause 
the wind to become proton rich. Under these conditions, 
the wind may synthesize some of th e light p-process ele- 
ments via the so called i^p-process ()Fr6hlich et al.|[200l 
iPruet et al.l[200l . 

After the initial deleptonization burst though, the 
net lepton number carried by the wind will be small. 
Since the anti-electron neutrinosphere sits deeper in the 
PNS than the electron neutrinosphere, the electron neu- 
trinos will be coole r than the electron anti-neutrinos 
(jWooslev et al.l[T994l ) . If this asymmetry is large enough, 
the wind can become neutron rich at later times. Com- 
bined with an a-rich freeze out, this can give conditio ns 
favorable for the r-process (IWooslev fc HoffmanlfT99l) . 

In both cases, the resulting nucleosynthesis is char- 
acterized by the integrated neutron to seed ratio after 
charged particle reactions freeze out. For proton-rich 
winds, alpha-particles recombine into ^^C by the stan- 
dard triple alpha reaction and then alpha-capture and 
proton-capture reactio ns carry the nuclear flow up to ap- 
proximately mass 60 (IWooslev fc Hoffmanl[l99l . The 
slowest reaction in this sequence is ''He(2a,7)"'^^C, so the 
total number of seed nuclei produced is equal to the num- 
ber of ^^C nuclei produced. When only free protons are 
present, the integrated neutron density is determined by 
the rate of anti-neutrino capture on free protons. An es- 
timate for the neutr on to seed ratio for proton-rich winds 
is given by equation IB4I A neutron to seed ratio of only 
a few is required to bypass the few long-lived waiting 
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Figure 1. Neutrino two-color plot produced using the analytic relations in the Appendix. A neutron star with gravitational mass 
1.4 Mq has been assumed with a neutrinosphere radius of 10 km. The total neutrino luminosity is assumed to scale as Lv^.tot = 
10^^ {(T^)/3.5 MeV)"* ergs"'^. This luminosity is split between neutrinos and anti-neutrinos so as to ensure that the net deleptonization rate 
of the PNS is zero. The thick black line corresponds to an electron fraction of Ye = 0.5. Above this line, neutron-rich conditions obtain 
and below it the matter is proton-rich. The white region is where there no free neutrons remain after charged particle reactions cease. The 
N = 50 (tan) region corresponds to final neutron-to-seed ratios between 0.01 and 15. The "first peak" (yellow) region corresponds to a 
neutron-to-seed ratio between 15 and 70, and the "second peak" (orange) region is where the neutron-to-seed ratio is greater than 70. The 
dashed lines correspond to the base ten logarithm of the mass loss rate in solar masses per second. 



oints which hind er production of the hght p-isotopes 
Pruet et all I2OQ60 . Still, it is challenging to produce 
even a small neutron to seed ratio, since the dynamical 
time scale of the wind is short compared to the anti- 
neutrino capture time scale. 

In the neutron-rich case, seed nuclei are produced 
by a different reaction sequ ence *He(an,7)^Be(Q!,n)^^C 
(jWooslev fc HoffmanI 119921) . For the conditions en- 
countered in the wind, the neutron catalyzed triple- 
alpha reaction proceeds about ten times as quickly as 
^He(2a,7)^^C. This increases the seed number compared 
with that obtained in a proton-rich wind with similar dy- 
namical properties. Also, since there are free neutrons, 
capture can proceed up to the N=50 closed shell iso- 
topes ^*Sr, *^Y, and ^°Zr. Here, the neutron density is 
just determined by the number of free neutrons left af- 
ter charged particle reactions freeze out. The neutron to 
seed ratio in neutron-rich winds can be approximated us- 
ing equation lBlOl Charged particle reactions continue up 



to the N=50 closed shell, at which point it becomes un- 
favorable to capture alpha particles due t o small separa- 
tion e nergies and large coulomb barriers (iHoffman et al.l 
|1996[) . If neutrons are exhausted during this process, the 
wind will mainly produce the isotopes ^ Sr, *^Y, and ^°Zr 
(iHoffman et al. 1997). This happens when the condition 



|. 0.42 -0.49.5^ 



1+1 



(1) 



is met. Here, /q w X^sj^a,! is the fraction of the initial 
helium abundance that gets processed into heavy nuclei. 
Neutron to seed ratios of approximately 30 and 110 are 
required to produce first and second peak r-process nu- 
cleosynthesis, respectively. 

Using the analytic results for the wind dynamics and 
nucleosynthesis given in the Appendix (equations IA6[ 
Ml rnil IA16IB11 [B4l [BToI and using the neutrino in- 
teraction rates given in ij3.1l to fix the thermodynamic 
state at Tc), one can easily explore the neutrino temper- 
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ature parameter space to determine the neutrino tem- 
peratures and fluxes that are most conducive to the r- 
process or the production of the light p-process. Figure 
[1] is a neutrino two-color plot where it is assumed that 
the deleptonization rate is zero and that the neutrino lu- 
minosity scales with the temperature to the fourth power 
{L^^^tot = 105i((Ti.)/3.5MeV)''ergs-i). The different 
nucleosynthetic regions are delineated by the final cal- 
culated neutron to seed ratio. To give a feeling for how a 
particular point in parameter space might contribute to 
the integrated nucleosynthesis of the wind, the mass loss 
rate is also shown. 

For a significant amount of material to move past the 
N = 50 closed shell during neutron-rich conditions, the 
anti-neutrino temperature must be approximately 60% 
higher than the neutrino temperature. For second peak 
r-process nucleosynthesis to occur, the asymmetry must 
be greater than 100%. Modern PNS co oling calculations 
do not give such larg e asymmetries (|Pons et al.l 119991 : 
iHuedepohl eral|[2009l i. 

Under proton-rich conditions, only a small region of the 
parameter space at high neutrino and low anti-neutrino 
temperature is favorable for the i^p-process. There will 
be a small amount of neutron production in the white 
region, but it is unlikely that significant production of 
the light p-process elements ^''Se, '^^Kr, ^"^Sr, and ^^Mo 
will occur. The region in neutrino temperature space 
where there is significant neutron production is unlikely 
to be reached. This region is small due to the short 
dynamical time scale of the wind, which reduces the time 
over which anti-neutrinos can capture on free neutrons. 
One should note that, very soon after shock formation in 
the supernova, a wind solution may not be appropriate 
and material will be entrained closer to the PNS for a 
longer period of time. Thi s scenario w o uld b e similar to 
the the conditions used in iPruet et"ari (|2006D . 

Therefore, based upon simple principles, it seems 
unlikely that the standard wind scenario will pro- 
duce r-process or light p-process isotopes in solar ra- 
tios, a s is required by ob servations of metal poor halo 
stars (ISneden et al.l[l99fih . This same conclusion has 
been reached by othe r authors ([Hoffman et al.l 119971 : 
[Thompson et al.l l200lD . but is repeated here in simple 
terms. We will find that our numerical calculations give 
similar results and that there is no significant r-process 
nucleosynthesis associated with the wind. Still, the wind 
can produce some isotopes that may have an observable 
signature. For standard PNS luminosities, the wind will 
spend a significant amount of time in the region of pa- 
rameter space were N = 50 closed shell nucleosynthesis 
occurs. 

3. COMPUTATIONAL METHOD 

To more accurately investigate the integrated nucle- 
osynthesis of the NDW, we have update d the implicit 
Lagrangia n hydrodynamic s code Kepler ([Weaver et al.l 
[1978; Wo oslev et al.|[200^ to carry out time-dependent 
simulations of the wind dynamics and nucleosynthe- 
sis. Kepler has b een used pre viously to study time- 
independent winds ([Oian fc Woos lev 1996), but the weak 
and nuclear physics employed there was rudimentary and 
nucleosynthesis was not tracked. Trajectories from Ke- 
pler were used for post-processing calculations of nucle- 
osynthesis in [Hoffman et al.. (,1997i ). 



Kepler solves the non-relativistic hydrodynamic equa- 
tions in Lagrangian coordinates assuming spherical 
symmetry. First order general relativistic correc- 
tions are included in the gravitational force law (cf. 
iSha piro fc Teukolskvl ([19831 )1 AU order v/c effects are 
neglected. This is justified since the maximum wind 
speeds encountered are, at most, a few percent of the 
speed of light. The momentum equation is then 



dvr 
~dt 



.dP Graf P 47rPr3\/ 2Gra\' 
dm r^ \ pc^ mc^ J \ rc^ J 

(2) 

where the symbols have their standard meaning s. As 
has be en shown by prev ious stud ies (lOian &: Woosleyl 
1996t ICardah fc Fulled [19971 : lOtsuki et al. l 120001 : 
Thompson et al.l I2001D . general relativistic correc- 



tions to the gravitational force can have an appreciable 
effect on the entropy and dynamical time scale of the 
wind. The equation of state includes a Boltzmann gas 
of nucleons and nuclei, an arbitrarily relativistic and 
degenerate ideal electron gas, and photons. 

3.1. Weak Interaction Physics 

Energy deposition from electron neutrino capture on 
nucleons, neutrino annihilation of all neutrino fiavors, 
and neutrino scattering of all flavors on electrons is in- 
cluded in the total neutrino heating rate. Neutrino 
"transport" is calculated in the light-bulb approxima- 
tion. The energy deposition rate is dominated by neu- 
trino capture s on nucleons . The neutrino annihilation 
rates given in iJankal ([19911 ) are employed. For the scat - 
tering rates, the rates given in lan fc WooslevI (IT996[ ) 
are used, but we include general relativistic corrections. 
Standard neutrino capture rates are employed in the 
limit of infinitely heavy nucleons with first order cor- 
rections. In this limit, the cross section is (Y.Z. Qian, 
private communication) 



(J„-n 



Gj,cos^{9c 
7r(ftc)4 



[5y + 3g^] (e. ± A)' (l±WM,^^e, 



(3) 

Here, Gp is the Fermi coupling constant, 9c is the 
Cabibo angle, gv and gA are the dimensionless vector 
and axial- vector coupling constants for nucleons, A is the 
proton neutron mass difference, is the neutrino energy, 
and Wm accounts for the weak magnetism and recoil 
corrections to the neutrino-nucleon cross section when 
the base cross s ection is derived in the limit of infinitely 
heavy nucleons ([Horowitdl2002| ) . This correction reduces 
the anti-neutrino cross section and increases the neutrino 
cross section (by about a total of 10% at the energies en- 
countered in NDWs), which, for a given incident neutrino 
spectrum, significantly increases the asymptotic electron 
fraction. Assuming a thermal distribution, these cross 
sections result in the neutrino energy deposition rate for 
anti-electron neutrino capture 



gpp = 4.2 X 10^**ergs s"^g 



-1 YpLp^^i 



-w: 



pp (4) 



(4) 



-(2A + W^|fA2)Q+A 
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and a similar expression for electron neutrino capture. 
The neutrino energy distributions are parameterized by 
assuming a Fermi-Dirac spectrum. The neutrino energy 
averages, (e"), are evaluated using this distribution. The 
neutrino energy moments and luminosity are evaluated 
in the rest frame of the fluid. With general relativistic 
corrections for the bending of null geodesies, the average 
neutrino angle is given by 

Special relativistic corrections are negligible in the re- 
gions where neutrino interactions are important. 

The lepton capture rates used are calculated in the 
limit of infinitely heavy nucleons. This results in a 
positron capture energy loss rate 

ge+„ = 6.9 X IQis ergs g"! s"! r„Tfo 

1 . 

here fdu, rf) — (exp(ii— r;) + l) , 77 is the electron degen- 
eracy parameter, 6 is the proton neutron mass difference 
divided by fc^T, and y„ is the neutron fraction. A similar 
rate is employed for electron capture. 

For the neutrino losses, we include electron and 
positron captur e on nucleons and include thermal losses 
as tabulated in lltoh et al.l ()1996[ ). The energy loss rate 
in the wind is dominated by the electron captures. 

3.2. Nuclear Physics 

During a hydrodynamic time step in Kepler, the nu- 
clear energy generation rate and the changing nuclear 
composition are calculated using a modified version o f 
the 19-isotope network described in lWeaver et al.l (|1978l ). 
Neutrino and electron capture rates on nucleons are cou- 
pled to the network, which are calculated under the 
same assumptions as the charged current energy de- 
position/loss rates described above. Therefore, non- 
equilibrium evolution of the electron fraction is accu- 
rately tracked. 

Although this network is appropriate for calculating 
energy generation throughout the entire wind, it is not 
large enough to accurately track the nucleosynthesis once 
alpha recombination begins at T « 0.5 MeV. Therefore, 
for temperatures below 20GK an adaptive network is run 
alongside the hydrodynamics cal culation. The details 
of this network can be fo und in iWooslev et all (|2004f ) 
and [Rauschcr et al.l (|2002l ) . As a fluid element passes 
the temperature threshold, the composition from the 19- 
isotope network is mapped into the adaptive network. 
Typically, the network contains approximately 2000 iso- 
topes. Where available, experimental nuclear reaction 
rates are employed, but the vast majority of the rates 
employed in th e network come from the statistical model 
calculations of iRauscher fc Thielemaiml (2000.' ) . In gen- 
eral, the nuclear physics employed in the s e calc ulations is 
the same as that used in IRauscher et al.l (|2002D . The nu- 
cleon weak interaction rates employed in the 19-isotope 
network are also used in the adaptive network. 

3.3. Problem Setup and Boundary Conditions 

To start the neutrino driven wind problem, an atmo- 
sphere of mass 0.01 Mq is allowed to relax to hydrostatic 



equilibrium on top of a fixed inner boundary at the neu- 
tron stars radius. The mass enclosed by the inner bound- 
ary is the neutron star's mass. The photon luminosity 
from the neutron star is assumed to be nearly Edding- 
ton, but we have found that the properties of the wind 
are insensitive to the the luminosity boundary condition. 
Once hydrostatic equilibrium is achieved, the neutrino 
flux is turned on and a thermal wind forms. This wind 
is allowed to relax to a quasi-steady state, and then the 
19 isotope network is turned on and the wind is, once 
again, allowed to reach a quasi-steady state. After this 
point, the neutrino flux is allowed to vary with time, and 
the adaptive network is turned on. 

As the calculation proceeds, the mass of the enve- 
lope being followed decreases and could eventually all 
be blown away. To prevent this, mass is added back to 
the innermost mass elements at a rate equal to the mass 
loss rate in the wind. The mass added to a fluid element 
at each time step is a small fraction of its total mass. We 
flnd that mass recycling has no effect on the properties 
of the wind. It is simply a way of treating a problem that 
is essentially Eulerian in a Lagrangian code. 

For most runs, a zero outer boundary pressure and 
temperature are assumed. To investigate the effect of a 
wind termination shock, a time dependent outer bound- 
ary condition is included in some of the simulations de- 
tailed below. The pressure of the radiation dominated re- 
gio n behind the superno va shock is approximately given 
by (IWooslev et al.l[200l 



where Esn is the explosion energy of the supernova, Vsn 
is the supernova shock velocity, and t is the time elapsed 
since the sho c k wa s launched. As was discussed in 
lArcones et al.l ()2007f ). this results in a wind termination 
shock at a radius where the condition Pwvti + Pw ~ Pps 
obtains, where is the wind velocity and pw is the 
wind density. To avoid an accumulation of too many 
zones, mass elements are removed from the calculation 
once they exceed a radius of 10, 000 km. This is well out- 
side the sonic point and nuclear burning has ceased by 
this radius in all calculations . 

4. NUMERICAL RESULTS 

To survey both low and intermediate mass core col- 
lapse supernovae, neutrino emission histories were taken 
from two core collap se calculations, one from a 2OAf0 
(jWooslev et al.l [Tool supernova calculatio n and the 
other from a 8.8Mq (|Huedepohl et al.ll2009l ) supernova 
calculation. Since the PNSs studied have significantly 
different masses and neutrino emission characteristics, 
one is able to get a rough picture of how integrated nu- 
cleosynthesis in the NDW varies with progenitor mass. 

4.1. Neutrino Driven Wind from a 20A/q Supernova 

The first set o f neutrino luminositie s and temperatures 
are taken from IWooslev et al.l (|1994) . This calculation 
began with a 20M (7) progenitor rneant t o model the pro- 
genitor of 1987A (jWooslev et"aIlll988D . The resuUing 
neutron star had a gravitational mass of 1.4Mq and the 
neutrino sphere was taken to be at 10 km. The neu- 
trino luminosities and average energies as a function of 
time from this model are shown in figure [2] After about 
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Figure 2. Neutrino luminosities and temperatures taken from 
the model of Wooslcy ct al. (1994). The top panel is the neutrino 
luminosities. The bottom panel is the average neutrino energies. 
The solid line corresponds to Ue, the dashed line corresponds to Ve, 
the dot-dashed line corresponds to v^^t- 

4 seconds, the neutrino energies become constant and 
the large difference between the electron neutrino and 
anti-neutrino energies implies that the wind will be neu- 
tron rich. This supernova model had some numerical 
deficiencies (Sam Dalhed, Private Communication). The 
entropy calculated for the wind in iWoosley et al.l ()1994f ) 
[S/NAk ~ 400) were unrealistically large due to some 
problems with the equation of state. Here, that is not 
so important because the NDW is being calculated sep- 
arately, but this study did rely on older neutrino inter- 
action rates and did not include weak magnetism cor- 
rections (see §3.1|1 . Therefore, the results obtained us- 
ing these neutrino histories are only suggestive of what 
might happen in a more massive star. If weak magnetism 
were taken into account, the calculated electron and anti- 
electron neutrino temperatures would probably be some- 
what further apart. 

The calculation was run for a total of 18 seconds. 
During this time, the mass loss rate decreased by al- 
most three orders of magnitude while a total mass of 
2 X 10~^ Mq was lost in the wind. A snapshot of the 
wind structure two seconds after bounce is shown in fig- 
ure[3] Note that the wind velocity stays very sub-luminal 
throughout the calculation. Therefore, the neglect of 
special relativistic effects is reasonable. The secondary 
bump in the energy deposition rate occurs at the same 
radius where nucleons and alpha-particles assemble into 
heavy nuclei. This increases the entropy by about 10 
units. Clearly, the electron fraction is set interior to were 
nuclei form. The radius where nu clei form is at a large 
enougn value that the alpha effect (jFuller fc Meverill995l ) 
is not significant at early times in the wind. However, as 
the neutrino luminosity decreases with time, nucleon re- 
combination occurs at a smaller radius, and the alpha 
effect becomes increasingly important. 

The time evolution of the wind as calculated by Ke- 
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Figure 3. Wind structure after two seconds in the model us- 
ing the neutrino luminosities from IWoosley et al.l l|1994l ). The top 
panel shows the density in units of 10** gcm~'^ (solid line) and the 
radial velocity in units of 10"^ kms~^ (dot-dashed line). The middle 
panel shows the net energy deposition rate from weak and strong 
interactions in units of 10 ergg~^s~^ (dot-dashed line) and the 
entropy (solid line). The bottom panel shows the temperature in 
units of 2 X 10^ K (solid line), the electron fraction (dot-dahsed 
line), and the fraction of material contained in nuclei (dotted line). 

pier is shown in figure SI The increase in asymptotic 
entropy is mainly driven by the decrease in neutrino lu- 
minosity, since the average neutrino energies do not vary 
greatly. The analytic approximation (calculated using 
equation IA8I and the neutrino interaction rates given in 
H3.1\i to the entropy tracks the entropy calculated in Ke- 
pler fairly well. This implies that the variation in the 
neutrino luminosity with time does not significantly al- 
ter the dynamics from a steady s tate wind. In contras t 
to the high entropies reported in iWooslev et al.l (|1994l) . 
the entropy here never exceeds 130. For the time scales 
and electron fractions also obtained, such a low value of 
entropy is not sufficient to give a strong r-process (see 
below). 

The electron neutrino and anti-neutrino energies do 
move further apart as a function of time though, which 
causes the wind to evolve from proton-rich conditions at 
early times to neutron-rich conditions later. A transition 
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Figure 4. Propert ies of the neutrino driven wind from tiie 
IWooslev et al.l 1)19941 ) supernova model as a function of time. Tiie 



thick Unes correspond to the numerical results from Kepler and the 
thin lines correspond to the predictions of the analytic estimates 
described in the appendix. The solid line is the dimensionless en- 
tropy per baryon, the dashed line is the electron fraction, the dash 
dotted line is the dynamical timescale, and the dotted line is the 
mass loss rate. All of the quantities are taken extracted from where 
the wind temperature reaches 2 GK. 



Figure 5. Isotopic productio n factors from the ND W model when 
the neutrino luminosities from IWooslev et al.l II1994I ) are used. The 
production factors are calculated assuming that 18.4 Mq of mate- 
rial was ejected in the supernova in addition to the wind. The top 
dashed line corresponds to the greatest production factor in the 
wind, the solid line is a factor of two below that, and the bottom 
dashed line is a factor of two below the solid line. These lines 
specify an approximate coproduction band for the wind alone. 



occurs from the synthesis of proton-rich isotopes via the 
lyp-piocess at early times to the a— process mediated by 
the reaction sequence a(Q!n,7)^Be(a,n)^^C later. The 
slight difference between the analytic approximation and 
the Kepler calculation of Yg is due to the alpha effect 
(jFuUer & Mover 1995). 

Integrated production factors for the wind are shown 
in figure El The production factor for the species i is 
defined as y i r 

where X^^u, is the mass fraction of species i in the wind 
after all material has decayed to stable isotopes, is 
the mass ejected in the wind, and Ms„ is the amount of 
mass ejected by the entire supernova. Xi^Q is the mass 
fraction of isot ope i in the sun for which the values of 
ILoddersI (fMol were used. The only isotopes that are 
co-produced in the wind alone are *^Rb, ^*Sr, *^Y, and 
^°Zr, with production factor of *^Sr about a factor of 3 
higher than the other two N = 50 closed shell isotopes. 
Before eight seconds, the production factors had been 
much closer. After eight seconds though, the wind is 
dominated by ^^Sr because Ye ~ 0.45 and only 53% of 
alpha particles are free after freeze out which puts ^ « 
0.41 of the heavy nuclei just below the range given in 
equation [TJ There are not enough free neutrons to make 
any significant amount of heavier nuclei, and this results 
in significant production of the stable N = 50 closed shell 
isotope with the lowest ^. 

During the first four seconds, the wind is proton rich 
and the isotopes ^SQa, ^^''^Ge, ^^Jegg^ 78,80,82j^j. ^^^^ 
produced by proton captures on seed nuclei produced by 
the triple-alpha reaction and subsequent {a,p) reactions. 
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Figure 6. Neutrino two-color plot when the anti-neutrino lumi- 
nosity is 1.2 times neutrino luminosity, and the total luminosity 
scales with average temperature to the fourth. Si milar to figure 
[T] Th e red line is the neutrino temperatures from IWooslev et al.l 
ll99lV 



Although the mass loss rate is much higher when the 
wind is proton rich, the alpha-fraction freezes out at 98% 
of its initial value, which results in significantly decreased 
production of heavy nuclei. The difference in final alpha 
fraction between the neutron- and proton-rich phases of 
the wind is due mainly to the difference in speed of the re- 
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Figure 7. Isotopic production factors from the NDW mod el em- 
ploying the neutrino luminosities from lWooslev et al.l 119941 ) with 
the anti-electron neutrino temperature reduced by 15%. The pro- 
duction factors are calculated assuming that 18.4 Mq of material 
was ejected in the supernova in addition to the wind. The horizon- 
tal lines are similar to those in figure [S] 
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Figure 8. Isotopic production factors from the NDW mod el em- 
ploying the neutrino luminosities from lWooslev et al.l l|1994l ) with 
weak magnetism corrections turned off. The production factors are 
calculated assuming that 18.4 Mq of material was ejected in the 
supernova in addition to the wind. The horizontal lines are similar 
to those in figure [5] 

action chains q;(2q!,7)^^C and a(an,7)^Be(a,n)^^C, but 
also to the decreased entropy at early times. 

We can compare this with the analytic predictions 
for nucleosynthesis by plotting the neutrino temperature 
evolution from this model on a neutrino "two-color plot" 
(figure [S]). Here we have set — 1.2L^ which is ap- 
proximately correct at late times in the calculation of 
IWooslev et al., (,1994 ). The wind never reaches a region in 



Figure 9. Isotopic productio n factors from the NDW model when 
the neutrino luminosities from lWooslev et al.l II1994I) are used and 
an external boundary pressure is specified as described in the text, 
which results in a wind termination shock. The production factors 
are calculated assuming that 18.4 Mq of material was ejected in 
the supernova in addition to the wind. The horizontal lines are 
similar to those in figure [5] 

which r-process nucleosynthesis is expected, but spends 
a significant amount of time making nuclei in the N = 
50 closed shell isotones. 

4.1.1. Variations in Neutrino Properties 

Since the neutrino temperatures from the original 
model were uncertain, several other models were calcu- 
lated. One had a reduced (by 15%) electron antineutrino 
temperature; another had the weak magnetism correc- 
tions to the neutrino interaction rates turned off. A 
smaller antineutrino temperature is more in line with 
recent cal culations of PNS cooling (iPons et al.l 119991: 
iKeil et all [2003. ). Because the model of IWooslev et al.l 
(1994) did not include weak magnetism corrections, our 
model with weak magnetism corrections turned off is 
more consistent with the original supernova model. 

The production factors for the model with a reduced 
electron antineutrino temperature are shown in figure [71 
The yield of ®®Sr is reduced by almost a factor of ten 
from the base case, while the production factors of 
and ^°Zr are reduced by a factor of three. In this case, 
the wind also produces the proton-rich isotopes ^*Se, 
^^Kr, and ^''Sr. The coproduction line for lighter ele- 
ments like oxygen in a 2OM0 supernova at solar metal- 
licity is around 18, so the wind could contribute to the 
total nucleosynthesis if the antineutrino temperature was 
reduced, but its contribution would be small. 

The yields when weak magnetism corrections are ig- 
nored are shown in figure |8l Without weak magnetism, 
the electron fraction drops below 0.4 at late times when 
the entropy is fairly high. Equation [1] is no longer sat- 
isfied and material moves past the N = 50 closed shell 
towards A ss 110. Some r-process isotopes are produced, 
such as ^^Zr and ^°°Mo, but not anywhere near solar 
ratios, and no material reaches the first r-process peak. 



Integrated Neutrino Driven Wind Nucleosynthesis 




Figure 10. Combined isotopic production factors of the neutrino dr iven wind with una l tered neutrino temperatures and including weak 
magnetism corrections added to those of a 20 Mq stellar model from lWooslev fc Weave3 1)19951 ). The solid black line is the coproduction 
line with ^®0. The dashed lines are a factor of two above and below the coproduction line. The neutrino driven wind is responsible for the 
production of **Sr, ^^Y, and ^OZr. 



4.1.2. Effect of a Wind Termination Shock 

To investigate the possible effect of a wind termination 
shock on nucleosynthesis, another model was run with a 
boundary pressure and temperature determined by equa- 
tion [71 An explosion energy of 10^^ erg was assumed and 
the shock velocity was taken as 2 x 10^ cm s""'^. This re- 
sulted in a wind termination shock that was always at a 
radius greater than 10'^ km. The production factors from 
this model are shown in figure |9l Similar to the simula- 
tion without a wind termination shock, the N=50 closed 
shell elements dominate the wind's nucleosynthesis. 

The main difference between the case with and without 
a wind termination shock is a shift in the mass of isotopes 
produced during the proton-rich phase. This can be seen 
in the increased production of Mo. During this phase, 
the post shock temperature varied from 2.5 GK down 
to 0.8 GK and the density varied from 5 x 10** g cm~'^ 
to 5 X 10^ g cm~^. These conditions are very favorable 
for continued proton capture once the long lived waiting 
point isotopes ^^Ni and ^^Ge are bypassed by (n,p) re- 
actions. Because these conditions persist for at least a 
second after a fluid element passes through the wind ter- 
mination shock, significantly more proton captures can 
occur on seed nuclei that have moved past mass ^ 64 
relative to the case with no termination shock. Still, not 
many more neutrons are produced per seed nucleus rela- 
tive to the base run. Therefore, the net number of seeds 
that get past the long lived waiting points remains small 
and the proton-rich wind does not contribute to the in- 
tegrated nucleosynthesis. It should also be noted that 
a different treatment of the wind's interaction with the 
supernova shock might result in a breeze solution which 
may supply m ore favorable c onditions for z/p-process nu- 
cleosynthesis (|Wanaidl2006[ ). 

4.1.3. Total Supernova Yields 

In figure [101 the pr oduction factors fr o m a 20M(7i 
supernova model from iWooslev fc Weaveil (|1995[ ) have 
been combined with the production factors we calcu- 
lated in the NDW wit h the unaltered neutrino histories of 
(jWooslev et al.|[T99^ with weak magnetism corrections 



included. The wind could be responsible for synthesizing 
the isotopes ^'^Rb, ^*Sr, *^Y, and ^"Zr. ^*Sr production is 
above the co-production band, but the rest are in agree- 
ment with the stellar yiel ds. This overproducti on of ^*Sr 
is similar to the result of iHoffman et al.l (|1997l ). 

For the model with a reduced anti-electron neutrino 
temperature combined with the yields from the 20Mq 
supernova model, the wind contributes 28%, 42%, 35%, 
75%, 75%, and 80% of the total ^^Se, ^^Kr, ^^Sr, s^Sr, 
*^Y, and ^"Zr abundances in the supernova, respectively. 
This wind model does not result in any isotopes being 
overproduced relative to the rest of the yields of the su- 
pernova. For the case with weak magnetism turned off, 
the nuclei produced by the wind are overproduced rel- 
ative to those made in the rest of the star by factor of 
nearly 100, hence this would need to be a very rare event 
if this model were realistic. 

Clearly, weak magnetism corrections and variations in 
the neutrino temperatures have a very significant effect 
on nucleosynthesis in the wind. Aside from the effects of 
an extra source of energy (|5.3p , the neutrino spectra are 
the largest current theoretical uncertainty in models of 
the NDW. 

4.2. Neutrino Driven Wind from a 8.&Mq Supernova 

The second PNS model is a more modern one- 
dimensional calc ulation of an electron-capture supernova 
(jHuedepohl et a l. 2009) that s tarted from an 8.8M0 pro- 
genitor model ( Nomotcilll984D . This resulted in a PNS 
with a gravitational mass of 1.27Mo and a radius of 15 
km. Together the lower mass and increased radius im- 
ply a lower gravitational potential at the neutrinosphere. 
This work employed neutrino interaction rates which 
took weak magnetism and "in-medium" effects into ac- 
count. The neutrino luminosities and average energies 
as a function of time are shown in figure 1111 The max- 
imum difference between the electron and anti-electron 
neutrino average energies is si gnificantly less than in the 
model of IWooslev et al.l (|1994[ ). This is likely due in part 
to both the decreased gravitational potential of the PNS 
and the more accurate neutrino interaction rates in the 
newer model. 
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Figure 11. Neutrin o luminosities and temperatures taken from 
tlie model of Hucd epohl et al.l II2009I ) . The line styles are the same 
as in figure [2] 

The calculation was run for a total of nine seconds, at 
which point the mass loss rate had dropped by two or- 
ders of magnitude. The total amount of mass ejected in 
the wind was 3.8 x 1O~*M0. In figure [T^ the proper- 
ties of the NDW calculated using Kepler are plotted as a 
function of time. Notice that the entropy never reaches 
above 100 in this model, which diminishes the likelihood 
of significant nucleosynthesis. For comparison, we also 
include the analytic estimates detailed above. There is 
reasonable agreement between the analytic and the nu- 
merical calculations, but not nearly as good as in the 
20Mq model. 

In contrast to the simulat i on run with the neutrino lu- 
minosities of iWooslev et al.l ()1994[ ). the electron fraction 
continues to increase with time. The difference between 
the average electron neutrino energy and electron anti- 
neutrino energy is, at most , about 3 MeV, comp ared to a 
maximum of 8 MeV in the lWooslev et al.l H1994D calcula- 
tions. Also, the difference between the average neutrino 
energies decreases as a function of time, com pared to 
an increase with time in IWoosley et al.l (|1994l ). Finally, 
the energies of a ll kinds of neutrinos are lower in the 
iHuedepohl et al.l 1^09) calculation, so that the proton- 
neutron rest mass difference significantly suppresses the 
anti-neutrino capture rate relative to the neutrino cap- 
ture rate. These differences are presumably due to both 
the different neutron star masses and neutrino interac- 
tion rates employed. 

The conditions in this model thus preclude any r- 
process nucleosynthesis, but they are potentially favor- 
able for production of some low mass p-process isotopes 
by the z^p-process. The integrated isotopic production 
factors are shown in figure[T3l The total ejected mass was 
take as 7.4 Mq, as 1.4 Mp, neutron s tar is left behind in 
the calculation of IHuedepohl et all (|2009[ ). During the 
calculation a maximum network size of 988 isotopes is 
reached. The p-process elements ^"^Se and ^*Kr are co- 



Figure 12. Pr operties of the neutrino driven wind from the 
IHuedepohl et al.l |2009;) supernova model as a function of time. 
The lines have the same meaning as in figure [4] 




Figure 13. Isotopic production factor s from the NDW model 
employing the neutrino luminosities from IH uedepohl et al.l 1)20091 ). 
The production factors are calculated assuming that 7.4 Mq of 
material was ejected in the supernova in addition to the wind. 
The horizontal lines are similar to those in figure \E\ Notice that 
none of the production factors are significantly greater than one. 

produced with ^^Cu, ^^Zn, and '^^Ga, but the maximum 
production factor for any isotope is 1 when weighted with 
the total mass ejected in the supernova. Therefore, in 
this simple model, the proton-rich wind from low mass 
neutron stars will not contribute significantly to galactic 
chemical evolution. 

The entropies encountered when the mass loss rate is 
high are low 50), so that there is more production of 
^^Ni by triple-alpha and a subsequent ap-process. As the 
neutron abundance available for the z/p-process is given 
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Figure 14. Isotopic production factor s from tlie NDW mode l 
employing the neutrino luminosities from lHuedepohl et al.l II2009I '). 
but including a time dependent external boundary pressure which 
results in a wind termination shock. The production factors are 
calculated assuming that 7.4 Mq of material was ejected in the 
supernova in addition to the wind. The horizontal lines are similar 
to those in figure [5] Notice that the production factors are almost 
unchanged when an external boundary pressure is added. 
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increased seed production reduces the available neutron 
abundance and therefore hinders production of the p- 
process elements ''''Se, ''^Kr, *'*Sr, and ^^Mo. Addition- 
ally, at early times, the dynamical time scale is short 
which implies a smaller integrated neutron to seed ratio, 
A„ (see the appendix). 

The yields of from this model cannot be combined with 
the yields from the r est of the supernova because they 
are not published. As iNomotol ([1984) has discussed, the 
mass inside the helium burning shell was close to the 
mass of the neutron star that was left after the explo- 
sion. Therefore the ejecta of the supernova is expected 
to have small production factors. This implies that, even 
when the yields of the NDW are combined with the rest 
of the supernova, it is unlikely that these low mass core 
collapse supernovae will contribute significantly to galac- 
tic chemical evolution. 

4.2.1. Effect of a Wind Termination Shock 

As was mentioned above, it is very possible that a tran- 
sonic wind solution may not be appropria t e this early 
in the supernovas evolution. iFischer et al.l (|2009[ ) have 
found that a wind termination shock is not present in 
a one -dimensi o nal su pernova model using the progenitor 
from ' Nomotol (1984). Still, it is interesting to consider 
the effect of a reverse shock on the wind nucleosynthesis. 

A second simulation was run with a time dependent 
boundary pressure given by equation [71 with Esn = 
2^q50 gj.g g^j^j = 2 X 10^ cm s^^. This results in a wind 
termination shock at a radius of approximately 3 x 10^ cm 
throughout the simulation. Inside the wind termination 



shock the wind dynamics are very similar to those in the 
run with no boundary pressure. The production factors 
from this model are shown in figure [141 Clearly, there 
is almost no difference in the nucleosynthesis in the runs 
with and without a wind termination shock. 

After 0.75 s, the post shock temperature drops below 
1 GK and the wind termination shock has little effect 
on subsequent nucleosynthesis. Because the post shock 
temperature is high for less than one second, the wind 
termination shock has very little effect on the integrated 
nucleosynthesis. A larger explosion energy would likely 
result in a larger effect on the nucleosynthesis, but there 
are still very few neutrons available to bypass the long 
lived waiting points and it seems unlikely that the pro- 
duction factors would be increased by more than a factor 
of a few. 

5. DISCUSSION 
5.1. Comparison with SN 1987 A 

Sinc e the progenitor model used in iWooslev et al.l 
()1994[) was a model for SN 1987A, it is interesting to 
compare our predicted abundances with those observed 
in the ejecta of that event. *^Sr, produced by the NDW, 
dominates the elemental strontium yield when the results 
of the 20 M(7) wind mode l are c ombined with those pre- 
dicted by IWooslev et al.l (|1988D , who ignored the wind. 
For the base NDW model, [Sr/Fe]= 0.8, if weak mag- 
netism corrections are neglected, [Sr/Fe]= 1.6; and if the 
anti-neutrino temperature is reduced in the base model 
b y 15%, fSr/Fe1= 0.2. 

iMazzaU etahl ()1992l ) found [Sr/Fe]?a 0.3 in the ejecta 
of SN 1987A 2-3 weeks after the explosion. This obser- 
vation has a substantial error bar due to the uncertainty 
of modelling the spectrum of the expanding ejecta. Of 
even greater concern is comparing of our models for bulk 
yields with supernova photospheric abundances observed 
on a given day when the observations are probably not 
even probing the innermost ej ecta. Still, if one ass umes 
that the [Sr/Fe] ratio found bv IMazzaU et al.l ()1992D rep- 
resents the value for all of the ejecta (for instance by as- 
suming that "mixing" was extremely efficient), a weak 
constraint can be put on the n eutrino fluxes and en - 
ergies predicted by the model of IWooslev et al.l (|1994l ). 
The NDW model with a reduced anti-neutrino tempera- 
ture is much closer to the observed value of [Sr/Fe] than 
the other two models. This suggests the anti-neutrino 
temperature may have been overestimated in the origi- 
nal model, a change that would be more consistent with 
more modern calculat ions of neutrino spectra formation 
in PNS atmospheres (jKeil et al.l 120031 ). But obviously, 
a meaningful constraint will require a more complete 
modeling of the multi-dimensional explosion and time- 
dependent spectrum of SN 1987A. 

5.2. Strontium and Yttrium in Halo Stars 

Since strontium and yttrium are abundantly produced 
in our models, it may be that the NDW has contributed 
to their production throughout cosmic history. An in- 
teresting possibility is that the abundances of these el- 
ements might trace the birth rate of neutron stars at 
an early time. Taking a standard r-process abundance 
pattern from meta l poor st ars with strong r-process en- 
hancments. iTravaglio et al.l (2004 ) find that 8% and 18% 
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of solar strontium and yttrium, respectively, are not pro- 
duced by either the "standard" r-process or any compo- 
nent of the s-process. It therefore seems plausible that 
charged particle reactions in the NDW could make up 
this "missing" component. 

Any nucleosynthesis that happens in the NDW will be 
primary, i.e. provided that the mass function of neu- 
tron stars at birth does not itself scale with metallicity, 
similar nucleosynthesis will occur for stars of any popula- 
tion. Below [Fe/H]~ —1.5, no component of the s-process 
contribut es to the abundances of N = 50 closed shell 
isotopes (jSerminato et al.l[2009h . If the NDW escapes 
the potential well of the PNS, and contributes to the 
galactic budget of N = 50 closed shell isotopes, it should 
provide a flo or to [Sr/Fe] and [Y/F e]. Based upon the 
arguments of iTravaglio et al.l (|2004l ). this floor would be 
at [Sr/Fe]w -0.18 and [Y/Fe]w -0.16. These numbers 
assume that when the main r-process source contributes 
in addition to the NDW, [Sr/Fe] and [Y/Fe] approach 
their solar values even though the s-process has yet to 
contribute. This is consistent with observations. 

In defining this floor, one must assume that the abun- 
dances in a particular star sample a large number of 
individual supernovae. This is because the production 
of N=50 closed shell elements likely depends on the 
PNS mass and therefore the progenitor mass. As we 
have found, [Sr/Fe]= 0.8 in the 20Mq model with re- 
duced anti-neutrino temperatures, but the 8.8M0 model 
produces no strontium. Observations show that be- 
low [Fe/H] 3, the spreads in [Sr/Fe] and [Y/Fe] in- 

crease significantly and the mean values falloff som e 
()Francois et al.|[2007l: iCohen et al.|[200l iLai et al.ll2008l ) . 
Single stars have values of [Sr/Fe] below the predicted 
floor. This could be because, at this metallicity, the 
metals in a particular star come from only a handful of 
supernovae. 

Another possible explanation of this variation is that 
supernova fall back varies with metallicity. Since the 
NDW is the innermost portion of the supernova ejecta, it 
will be the most susceptible to fallback. It has been found 
that the amount of supernova fallback depends strongly 
on the metallicity of the progen itor, especiall y going be- 
tween zero and low metallicity (jZhang et al.l 120081. Ad- 
ditionally, mixing is also greatly reduced in zero metal- 
licity stars compared to sol ar metallicity stars du e to the 
formers compact structure (jJoggerst et al.ll2009| ). 

The current understanding of supernova fallback sug- 
gests that the nucleosynthetic contribution of the NDW 
will be suppressed at very low metallicity. Of course, the 
ejection of iron by the supernova is also very suscepti- 
ble to fallback, so the effect of fallback on the evolution 
of [Sr,Y/Fe] is complicated and may require flne tuning 
to give the observed dec rease. A somewhat d ifferent ex- 
planation was offered bv lOian fc Wasserburd ()20G8.) who 
attributed the fall off of [Sr/Fe] at low metallicity to the 
evolution of the "hypernova" rate with metallicity. For 
their purposes, hypernovae were stellar explosions that 
contributed iron without making much strontium. 

Given the sensitivity of strontium and yttrium yields 
to uncertain NDW characteristics, especially neutrino 
fluxes and temperatures, it may be some time before the 
complex history of these elements is even qualitatively 
understood. It is likely though that their abundances in 
halo stars will ultimately be powerful constraints upon 



the evolution of supernova physics as a function of metal- 
licity. 

5.3. Possible Modifications of the Basic Model 

As is clear from flgures [S] and [T31 the simplest case of 
a non-magnetic non-rotating NDW from a neutron star 
without additional energy deposition does not produce 
r-process nuclei in significant abundances. Are there ex- 
tensions to this simple scenario that could make the wind 
a site of the r-process? 

As was pointed out by IMetzger et al.l (j2007f) , the com- 
bination of rotation and magnetic fields can decrease the 
dynamical time scale by magnetic "flinging" . This is not 
particularly effective. Adding a non-thermal source of ki- 
netic energy means that less thermal energy must be put 
into the wind for it to escape the potential well. There- 
fore, lower entropies are achieved. It seems unlikely that 
this mechanism, by itself, will salvage the NDW as a site 
for the full r-process. If there were a way to make the 
rotation rate of the PNS high enough, it might be pos- 
sible that there would be a centrifugally driven outflow. 
Then the electron fraction would be determined by ki- 
netic equilibrium much deeper in the PNS envelope, and 
the material in the outflow would have an electron frac- 
tion much lower than that seen in the wind. 

To test this possibility, we ran calculations with a cen- 
trifugal force term added and corotation with the PNS 
enforced out to 10^ km. Unfortunately, for reasonable 
PNS spin rates (20 ms period), we found this had little 
effect on the nucleosynthesis. These calculations were 
in a regime were the electron fraction was still set by 
neutrino interactions. 

Many authors have disc ussed the possible ef- 
fects of both matter-enhanced(!Oian & Fuller' 1995'; Sid 
1995f) and collec tive neutrino (Pastor & Raffclt 2002( 
Duan et"aril2006t ) oscillations on NDW nucleosynthesis. 
If electron antineutrinos could undergo a collective oscil- 
lation near the launch radius while the electron neutrinos 
did not, this would increase the average energy of the an- 
tineutrinos if the fi and r neutrinos have a significantly 
higher temperature, facilitating a reduction in the elec- 
tron fraction. For a normal mass hierarchy however, mat- 
ter enhanced neutrino oscillations would probably cause 
electron neutrino flavor conversion, which would increase 
the electron fraction and decrease the probability of si g- 
niflcant r-process nucleosynthesis (|Oian fc Fulleij lT995^. 

Collective neutrino oscillations can cause antineutrino 
oscillations in the region were the electron fraction is set, 
and thereby decrease the electron fraction where pure 
MSW oscillat ions would have p redicted an increased elec- 
tron fraction( |Duan eFaDEOOi). Clearly, the main effect 
of oscillations would be on the composition of the wind, 
not the dynamics. As can be seen in the neutrino two 
color plots, oscillations would have to change the effective 
temperature of the anti-neutrinos by a very large amount 
to move from a region where N=50 close shell nucleosyn- 
thesis occurs to a region where the second r-process peak 
can be produced. 

These effects are based upon the assumption that /x- 
and T-neutrinos are signiflcantly more ene rgetic than the 
electron neutrinos. In the calculation of iWooslev et al.l 
(|1994D . this is the case, as can be seen in figured] In- 
terestingly, the fj, and r temperatures are almost the 
same as the electron anti-neutrino temperature in the 
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iHuedepohl et alj (|2009D calculation, which can be seen 
in figure [11] It is not clear wether this difference ob- 
tains because of the difference in the PNS masses or the 
significantly different neutrino physics employed in the 
calculations. A detailed study of neutrino transport in 
static backgrounds showed that the inclusion of all rele- 
vant neutrino interactions brings the average energies of 
the IX- and r- neutrinos closer to the tem perature of the 
anti-electron neutrinos ()Keil et al.ll2003() . Therefore it is 
uncertain wether or not neutrino oscillations could effect 
nucleosynthesis significantly. Clearly, the uncertainties 
here are not in the wind itself but in the formation of 
the spectra in the PNS and the details of neutrino trans- 
port with neutrino oscillations. 

Finally, i t has been suggested (Qian fc WooslevI [19961: 
ISuzuki fcN agataki 2005) that adding a secondary source 
of volumetric energy deposition can significantly increase 
the entropy of the wind, which results in a more alpha- 
rich freeze out and conditions that would be more favor- 
able for r-process nucleosynthesis. The addition of en- 
ergy to the wind also decreases the dynamical timescale. 
Since the important quantity to consider for the r-process 
is /tcl ((Hoffman ct al. 1997), both effects increase the 
chance of having a significant neutron to seed ratio after 
freeze out. If the NDW model is to be salvaged, this 
seems to us the minimal necessary extension. Of course, 
the physical process contributing this extra energy is very 
uncertain. One possibility is that oscillations of the PNS 
power sound waves which produce shocks and deposit 
en ergy in the wind, sirn ilar to the supernova mechanism 
of [Burrows et all ([2006[ ), but smaller in magnitude. We 
will explore this possibility in some detain in a subse- 
quent paper. 

From a chemical evolution standpoint, it is impor- 
tant to consider what effect neutron star mergers will 
have on the evolution of the r-process abundances. It 
seems unavoidable that r-process nuclei will be pro- 
duced in the tidal tails ejected during these mergers 
([Freiburghaus et "all [19990 and the amount of material 
ejected in these events is approximately enough to ac- 
count for the galactic inventory of r-process elements 
given the expected merge r rate ([Lattimer fc Schramm! 
[I976t IRosswog et al.l[l999l) . 

Neutron star mergers have been largely discounted be- 
cause inferred merger rates are small at low metallicity 
due to the long in spiral time and therefore they cannot 
account for the r -process enrichmen t seen in low metal- 
hcity halo stars ([Argast et al.ll2004H . Of course, the in- 
ferred merger rates are very uncertain as are models for 
the early evolution of the Milky Way, so that both the 
conclusion that neutron star mergers can account for the 
r-process inventory of the galaxy and that they are not 
consistent with producing the r-process at low metallic- 
ity are very uncertain. Clearly, there is significant room 
for more work in this area. 

Therefore, it seems reasonable that the galactic r- 
process abundances could be accounted for by a combina- 
tion of mergers and winds with an extra source of energy, 
either from an acoustically and/or magnetically active 
PNS. Since not every supernova will have the requisite 
conditions for an r-process, there will be significant vari- 
ation in the yields from single supernovae. This, along 
with the contribution from neutron star mergers, will 
give significant variation in the [r-process/a-element] val- 



ues found in single stars at low metallicity but averaged 
over many stars these should track one another, w hich is 
consistent with observations ([Sneden et al.|[2008[ ). 

6. CONCLUSIONS 

We have performed calculations of the dynamics and 
nucleosynthesis in time dependent neutrino driven winds. 
This was done for two sets of neutrino spectra calculated 
in one-dimensional supernova models taken from the lit- 
erature. The nucleosynthesis in these models was com- 
pared with supernova yields to determine if these mod- 
els were consistent with observations. Additionally, we 
compared the results of these numerical models to ana- 
lytic models of the neutrino driven wind and found good 
agreement. 

Similar to most o f the work on the NDW after 
IWooslev et all ([Tool . we find that it is unlikely that 
the r-process occurs in the neutrino driven wind unless 
there is something that causes significant deviation from 
a purely neutrino driven wind. Additionally, in the sim- 
plest case, there is little production of p-process elements 
at early times in the wind. In our calculation that used 
spectra from a more massive neutron star, the wind only 
produces the N=50 closed shell elements ^^Rb, *^Sr, ^^Y, 
and s^Zr. 

This result is sensitive to small changes in the neu- 
trino interaction rates (i.e. the inclusion of weak mag- 
netism) and changes to the neutrino temperature of order 
10%. Comparing our models with the over abundance of 
strontium seen in SN 1987A suggests that the difference 
between the electron a nd anti-electron neutrino temper- 
atures in the model of [Wooslev et al.[ ([19941) may have 
been to large. We also find that the effect of a wind 
termination shock on the wind nucleosynthesis is small. 

Using neutrino spectra from an S.SMq supernova that 
drive s a wind which is proto n rich throughout its dura- 
tion C Huedepohl eFall [20091) . we find that no significant 
i^p-process occurs and the wind docs not contribute to 
the yields of the supernova. The neutrino spectra from 
this model are pro bably more accurate than the spectra 
from the model of IWooslev et al.l ([1994D . We also inves- 
tigated the effect of an outer boundary pressure which 
resulted in a wind termination shock. This had a negli- 
gible effect on the nucleosynthesis. 

However, one also expects that the nucleosynthesis in 
the NDW will vary considerably from event to event, es- 
pecially with the mass and possibly the rotation rate of 
the PNS. The winds from more massive PNS have greater 
entropy and might, in general, be expected to produce 
heavier elements and more of them. The neutrino spec- 
tral histories of PNS as a function of mass have yet to be 
determined over a wide range of parameter space. Cur- 
rently, the neutrino luminosities and temperatures are 
the largest uncertainties in models of the NDW. 
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APPENDIX 
ANALYTIC WIND DYNAMICS 



To understand the wind dynamics, we follow arguments similar to those of lOian & WoosleV (19M) and 
ICardall &: Fuller, ((1997 ). Conservation of the stress-energy tensor and number flux in a Schwarzschild geometry in 
steady state leads to the wind equations in critical form 



M — Anr rribnjyv 
d . , , , quih 



jyv—\ii{-fyhr) - ^2 



C^hr 
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here M is the rest mass loss rate, mi, is the baryon mass, n is the baCTgn number density, v is the velocity measured 
by an observer at rest in the Schwarzschild frame, 7 = (l — w^/c^) , q is the neutrino heating rate per mass, Cg 
is the sound speed, M is the neutron star mass, hr — I + e/{mbC^) + P/{nmi,c^), P is the pressure, and e is the 
energy per baryon not including the rest mass. To fully describe the wind, these equations must be supplemented by 
a set of nuclear rate equations and an equation of state. For our analytic calculations we used a radiation-dominated 
non-degenerate equation of state comprised of relativistic electrons, positrons, and photons. 

First, we estimate where the critical radius (i.e. the radius where radiation pressure equals the nucleon gas pressure) 
sits in relation to the neutrino sphere. Neglecting the temperature gradient in the equation of hydrostatic balance 
gives the density structure of the atmosphere, 

GM mt, 



log(n(r)/nc) 



The optical depth of this atmosphere for neutrinos is 
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Taking the neutrinosphere to be at an optical depth of 2/3, the gas pressure equal to the radiation pressure at the 
critical radius, and approximating gravity as constant throughout the envelope, we arrive at an equation for the critical 
radius 



rc~Rv 



grubRt, \3to 



(A6) 



where tq = av{e^)ncTc/ gmb- For characteristic values of Ljy, e^, R^, and M, rc is only a few percent larger than 
R^. This implies that the GR corrections to the neutrino interaction rates at r^ will be at most a few percent. For 
characteristic values, the GR correction to gravity will be 2/(n)~^ ~ y^^ ~ 1-5. This agrees with the observation 
made by previous authors that GR correction s to gravity dominate over corrections to the neutrino interaction rates 
(ICardall fc Fulled [19971 iThompson et al.ll200lD . 

Assuming that most neutrino heating occurs near r^, the entropy can be considered constant once the temperature 
cools to the nucleon recombination temperature (T « 0.5 MeV). Therefore, the final nucl ear abundances in the 
wind depend mainly on the wind entropy, electron fraction, and the timescale for outflow (iQian fc Wooslevl[l99fil : 
iHoffman et al.llT997t) . To determine the contribution of the wind to the nucleosynthesis of the entire supernova, the 
mass loss rate must also be known. We now find estimates for these quantities and for the transonic radius of the 
wind. 

To determine the asymptotic entropy, the total energy deposition per baryon needs to be estimated. Using equation 
~]and assuming the asymptotic velocity is small, the total energy deposited per baryon is 

\n{jf) ~ \n{y,h,) « - \n{y,h,) « / dr-^^ « / dr-^^^ = Q/{mbC^) (A7) 
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Considering that most of the neutrino energy is deposited near the hydrostatic atmosphere, the final entropy per 
baryon is approximately 
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Assuming that the neutron-proton rest mass difference is neghgible, the entropy of the envelope is negligible, and 
taking the relativistic enthalpy outside the logarithm to be one results in the scaling relation 

Sf « 464 In {y-') ^'i.li'^.l^eV (|) ■ (A9) 

Notice that all the general relativistic corre ctions, both due to gravity and to the neu trino interaction rates, increase 
the entropy from the non-relativistic case (ICardah fc Fulled [l997l: lOtsuki et al.l[2000i) . As was discussed above, the 
dominant correction to the newtonian case is from the GR correction to gravity. 

To fix the mass loss rate, we must estimate the velocity at the critical radius. Taking the momenturn equation in 
critical form and assuming approximate hydrostatic equilibrium and subsonic velocities gives c^dv/dr w Assuming 
that acceleration has occurred over a scale height, we have 

~ ^ (AlO) 

where the heating rate at the critical radius is divided by two to account for the fact that beneath the critical radius 
the net heating goes to zero over approximately a scale height so that a characteristic value of the heating rate is one 
half the heatii ig ra te at the critical radius. The scale height is given by h~jj: — (p + P)GM/Pr^y^ . Combining these 
with equation lAll results in a mass loss rate of 

M«47rr4-^y2!W^ 



Using our result for the entropy gives a scaling relation for the mass loss rate, 

r4.i5/?.fl"/3 ., /„, \ 10/3 

Rlf{M/M^)ln{yi'y' W 

Notice that GR corrections reduce the mass loss rate significantly, by about a factor of 2. Only including GR effects 
on neutrino propagation and energies decreases the mass loss rate by about 10%. The reduced mass loss rate due to 
GR corrections does effect the integrated nucleosynthetic yields, although the effect is not as great as the effect of the 
GR corrections on the entropy. 

Given the entropy and the mass loss rate, we can solve for the evolution of temperature with radius outside of the 
heating region where the velocity is still small using the relation ln(/iry) — ln(2/c^r,c) + Q ~ 0, which gives 



2GM 



I - {l+Tsf/niby^ (A13) 



The dynamical timescale of the wind is defined by — vjy /n\dn/dr\. Steady state baryon number conservation 
yields 



dn 
dr 



2 Y dv CM 
r V dr ir-2y2 
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Neutrino energy deposition will no longer dominate the momentum equation when nucleons reform, but the velocity 
will be subsonic. In this limit, the momentum equation yields 

V dr r-^y^c^ r 
which results in an estimate for the dynamical timescale far inside the sonic point, 

GMM 
Awr^nmby^c^ 

For seed formation we are interested in the dynamical timescale around 2 MeV. Combining with our result for the 
temperature structure of the atmosphere we have 



16TT{GMmbymbK ^ 

r.a^ y (^17) 



To agree with the definition of the dynamical timescale given in lQian fc WooslevI (Il996f). this should be multiplied 
by three as our definition of differs slightly from the one used in IQi an fc WooslevI (I1996D . We note that all of the 
scaling relations above are equivalent to those of iQian fc Wooslevi (|1996l ) in the non-relativistic limit. 
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From this discussion, it is unclear if a transsonic wind will obtain. A reasonable criteria for transsonic solutions is 
that pv^ + P at the sonic poi nt is greater than the pressure behind the supernova shock, which is approximately given 
by equation [T) Equation IA3I along with equation of state for the wind can be combined to give the temperature and 
sound speed at the sonic radius in terms of the sonic radius and known quantities 



This results in an implicit equation for the sonic radius 



.V ~ ,-3/2 M""^'^' .. , [ 'iGM„smt I 



7/2 



where K is the radiation constant for photons and leptons combined. In the non-relativistic limit, the sonic radius 
reduces to 

f Sf \-»/^ / Mns V^^ f M \^^'^ , 
860km (^) (iFHMe^J ■ ^^'^^ 

The wind termination shock position is approximately given by ()Arcones et al.|[2"007l ) 

where Esn is the supernova explosion energy, Rs is the radius of the supernova shock, and is the wind velocity just 
inside the wind termination shock. 

ANALYTIC WIND NUCLEOSYNTHESIS 

The electron fraction in the wind is given by kinetic equilibrium of neutrino interaction rates at th e critical radius, 
as the temperature has decreased enough that lepton capture is unimportant (jQian fc Wooslevl[T996l) 

where X^^ is the neutrino capture rate per baryon and Ap_, is the anti-neutrino capture rate per baryon. After weak 
interactions cease and the temperature has decreased about 0.5 MeV, alpha particles form in the wind. The initial 
alpha number fraction is Y^^i ~ for neutron-rich conditions and Y^^i w 1/2 — Yf^/2 for proton-rich conditions. 

In both proton and neutron-rich winds, the nucleosynthesis will be characterized by the neutron to seed ratio in the 
wind. For winds with Y^ >« 0.5, alpha particles recombine into ^ ^C by the standard triple alpha reaction and then 
experience alpha particle captures up to approximately mass 56 (jWooslev fc Hoffmanl 11992!). The slowest reaction 
in this sequence is "^116(20,7)^^0, so the total number of seed nuclei produced is equal to the number of ^^0 nuclei 
produced. The rate of alpha destruction is given by 

^ « -UplY^h. (B2) 

where po ~ mi,n is the rest mass density and X^a is the rate of triple alpha, which includes double counting factors. 
The factor of 14 comes from assuming alpha captures stop at '^^Ni. In general, we define an abundance by rii = nYi, 
where rii is the number density of species i and n is the baryon density. Using our definition of the dynamical timescale, 
we have dr « —^dT. Transform ing Yg to a function of temperature makes the integral given above trivial to solve. 
Using a rate for triple-alpha from lOaughlan fc Fowled ()1988() and assuming the reaction fiow stops at ^^Ni results in 
a seed abundance at the end of the a-process given by 



28 



1 + 1.4X 10Vds72(i_y^)2l (B3) 



where sj^ enters because po / Sf for a radiation dominated equation of state and the density enters to the second 
power. 

Under proton-rich conditions, the j/p-process has the potential to occur. This process is similar to the rp-process, 
except that long lived beta-decays are bypassed by (n,p) reactions. An estimate for the integrated number of free 
neutrons produced is TdAp^(Tg w 2), so that the neutron to seed ratio is 

A« -TdA^.lTg «2)|^. (B4) 

Here, Ap. (T9 w 2) is the neutrino capture rate at the seed formation radius. A similar relation is found in lPruet et al.l 
(|2006D . Although it is hard to estimate its effect, we note that the presence of a reverse shock can significantly affect 
the z^p-process nucleosynthesis, as passage through the reverse shock slows the outward flow and rarefaction of the 
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wind. Additionally, it increases the temperature to close to the post supernova shock temperature. At early times for 
characteristic explosion energies, the wind is shock heated to a temperature of a few GK. This all combines to give a 
longer period of time over which proton capture on heavy nuclei is efficient and allows the z/p-process to continue to 
higher mass than it would if no wind termination shock were present. 

In the neutron-rich case, seed nuclei are produced by the slightly different reaction sequence ^IIe(an,7)^Be(a,n)^^C 
IWooslcy & Hoffman 1992). For the conditions encountered in the wind, neutron catalyzed triple-alpha proceeds 
about ten times as quickly as ^He(2a,7)^^C. This implies that there will be a larger seed number than in proton-rich 
conditions. The rate of helium destruction is given by the equations (jHoffman et al.|[l997[ ) 

rIY 7 

- -^P^Y^Y,Xo.,,CBe) (B5) 
ar 2 

^^ = 32^n^j >^nrip2T-3exp((i?9-2B„)/T) (B7) 

where Z is the average proton number of the seed nuclei, A is the average nucleon number of the seed nuclei, Na 
is Avogadro's number, is the alpha particle abundance, y„ is the neutron abundance, Ig is the ^Be abundance, 
Aa,n(^Be) is the rate of ^Be destruction by (a, n), and and i^g are the binding energies of ^He and ^Be, respectively. 
For Aa,n(^Be), we employ the rate given in lWrean et al.l ()1994[ ). This set of equations can be solved analytically by once 
again transforming from proper time to temperature using the dynamical timescale. The resulting implicit expression 
for the final alpha fraction is somewhat cumbersome, so we do not reproduce it here. The seed abundance in terms of 
the initial and final alpha fraction is v _ V 

Y, = 2 - "'^ (B8) 

the final neutron fraction is given by 



so that the neutron to seed ratio is 



A — 27 

Ynj = (1 - 2Y,) - 2^— - Y^j) (B9) 

The seed abundance at the end of charged particle reactions can be estimated by (jHoffman et al.|[l997f) 

Ys « (l - exp [-8 X lO'r.sfY^'] ) (Bll) 

where it has been assumed that the neutron abundance is what limits the reaction. Notice that the seed abundance 
in the proton-rich case depends on the entropy squared because it is mediated by an effective three body reaction, 
but in the neutron-rich case the entropy enters to the third power because of the effective four body interaction that 
mediates seed production. 
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